NASA TM X-182 


NASA TM X-182 





TECHNICAL MEMORANDUM 

X-182 


SOME EFFECTS OF VARIATIONS IN DENSITY AND AERODYNAMIC 

PARAMETERS ON THE CALCULATED FLUTTER CHARACTERISTICS 

OF FINITE-SPAN SWEPT AND UNSWEPT WINGS AT 

SUBSONIC AND SUPERSONIC SPEEDS 

ByE. Carson Yates, Jr. 

Langley Research Center 
Langley Field, Va, 


Declassified February 6, 1962 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON January 1960 




S=- CTs .pr t ” 1 


IX 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
TECHNICAL MEMORANDUM X-182 

SOME EFFECTS OF VARIATIONS IN DENSITY AND AERODYNAMIC 
PARAMETERS ON THE CALCULATED FLUTTER CHARACTERISTICS 
OF FINITE -SPAN SWEPT AND UNSWEPT WINGS AT 
SUBSONIC AND SUPERSONIC SPEEDS 
By E. Carson Yates, Jr. 

SUMMARY 


Subsonic and supersonic flutter calculations have been made by the 
method of NACA Research Memorandum L 57 LIO for 12 swept and unswept wings 
in order to evaluate the effects on calculated flutter characteristics of 
variations in aerodynamic parameters and flow density. Results are com- 
pared with experimental flutter data. With linearized-theory aerodynamic 
parameters used as the basis, the variations of aerodynamic parameters 
consisted of systematic constant-percent-chord changes in local aerody- 
namic center and systematic changes in the magnitude (but not in the 
distribution) of section lift-curve slope. These calculations indicated 
that when local aerodynamic centers are close to local centers of gravity, 
flutter calculations become extremely sensitive to small changes in 
aerodynamic -center position and calculations based on linearized- 
aerodynamic theory are likely to be considerably in error. If the local 
aerodynamic centers move rearward of the elastic axis, it is possible to 
have flutter speed increasing with increasing lift-curve slope. When 
local aerodynamic centers are sufficiently forward of the elastic axis 
and of the local centers of gravity, calculated flutter speed for a given 
wing varies approximately as the inverse square root of the lift-curve 
slope and as the inverse square root of the distance between aerodynamic 
center and center of gravity. 

Flutter calculations made for the 12 wings at several values of flow 
density indicated that flutter at constant dynamic pressure was indicated 
only for unswept wings for which local aerodynamic centers were not close 
to local centers of gravity. 
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INTRODUCTION 


Reference 1 presented a strip-theory type of flutter calculation 
procedure for finite-span swept and unswept wings based on spanwise dis- 
tributions of lift and pitching moment derived from distributions of 
aerodynamic parameters associated with the undeformed wing in steady 
flow . 1 In references 1 and 2 comparisons of experimental flutter data 
for several wings with flutter characteristics calculated for these wings 
by the modal -analysis method of reference 1 indicated that the method 
yielded accurate flutter results for a broad range of wings at subsonic, 
transonic, and supersonic speeds. However, for some of the unswept wings 
of reference 1, calculated flutter speeds were high compared with those 
obtained by experiment at supersonic Mach numbers. Since these calcula- 
tions utilized steady-flow aerodynamic parameters obtained from linearized- 
aerodynamic theory, it was hypothesized that for these unswept wings the 
absence of close agreement between experimental flutter speeds and those 
calculated from linearized-theory aerodynamics was related to the close 
proximity of the local aerodynamic centers to the local centers of gravity 
at supersonic speeds and to the fact that linear theory predicts too-far- 
rearward aerodynamic centers. This hypothesis was supported in reference 2 
by the greatly improved agreement between experimental flutter character- 
istics and those calculated by using steady-flow aerodynamic parameters 
obtained from flight test. (See figs. 16 and 18 of ref. 2.) The local 
aerodynamic centers obtained from the flight-test load distributions were 
characteristically forward of those predicted by linear theory. 

In view of this demonstrated inadequacy of linear! zed -aerodynamic 
theory and in view of the apparent sensitivity of some flutter calcula- 
tions to small changes in aerodynamic parameters, it was considered desir- 
able to obtain some quantitative indications of the effect on calculated 
flutter characteristics of variations in the aerodynamic parameters used. 
Accordingly, flutter calculations are presented herein for some of the 
swept and unswept wings of references 1 and 2 in which systematic modifi- 
cations have been applied to the linearized-theory aerodynamic parameters 
of reference 1. Specifically, two sets of subsonic and supersonic flutter 
calculations have been made for three wings of aspect ratio 4.0, taper 
ratio 0.6, quarter-chord sweep angle 45 ° > and varying chordwise center-of- 
gravity positions and for four unswept wings of aspect ratios 4.0 and 7*^ 

4n the method of reference 1, spanwise distributions of steady-flow 
section lift -curve slope and local aerodynamic center for the undeformed 
wing are used in conjunction with the "effective" angle -of -at tack distri- 
bution resulting from the assumed vibration modes in order to obtain values 
of section lift and pitching moment. Circulation functions modified on 
the basis of loadings for two-dimensional airfoils oscillating in com- 
pressible flow are employed to account for the effects of oscillatory 
motion on the magnitudes and phase angles of the lift and moment vectors. 
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and taper ratios of 0.6 and 1.0. The first set of calculations employed 
spanwise distributions of section lift-curve slope obtained from linear 
theory and spanwise distributions of linear-theory local aerodynamic center 
to which various constant-percent-chord increments had been added. The 
second set of calculations employed distributions of linear-theory section 
lift -curve slope which had been multiplied by various constant factors and 
distributions of linear-theory local aerodynamic centers. The first set 
of calculations thus indicates the sensitivity of calculated flutter char- 
acteristics to arbitrary constant -percent-chord changes in local aero- 
dynamic center, whereas the second set indicates sensitivity to arbitrary 
changes in the magnitude (but not in the distribution) of local lift- 
curve slope. 

The magnitude of section lift-curve slope and the position of local 
aerodynamic center are, of course, not arbitrarily variable on a physical 
wing. Indeed, these parameters are uniquely determined by wing geometry 
and flow conditions. The two sets of calculations mentioned previously, 
therefore, do not necessarily represent physically attainable conditions. 
They do, however, give some indication of the reliability of linear-theory 
results and may also be used to obtain an indication of possible effects 
on flutter characteristics of finite-thickness airfoil shapes or other 
factors which might change aerodynamic parameters from the linear-theory 
values. 

Subsonic and transonic flutter experiments (refs. 3 to 5) have indi- 
cated that flow density affects wing flutter characteristics in such a 
way that for a given wing at a given Mach number flutter generally occurs 
at approximately constant dynamic pressure, at least over a range of mass 
ratio (ratio of fluid mass to wing mass) normally encountered in flight 
and wind-tunnel testing. However, because of the number of variables 
involved, it is difficult to determine accurately by experiment the varia- 
tion of flutter characteristics with density. The effects of density 
variation on flutter characteristics have previously been calculated by 
using two-dimensional incompressible-flow aerodynamic parameters. (See 
refs. 3 to 3; for example.) As indicated in the discussion of reference 
these calculations may be expected to indicate flutter at approximately 
constant dynamic pressure only for high-aspect-ratio swept wings or for 
unswept wings. Inasmuch as the two-dimensional incompressible-flow type 
of flutter calculation used in reference 5 can be considered as a special 
case of the flutter calculation method of reference 1, it is considered 
desirable to examine the effect of density variations on flutter charac- 
teristics calculated by the method of reference 1 for several finite- 
span swept and unswept wings at subsonic and supersonic Mach numbers. 
Accordingly, flutter characteristics calculated by the method of refer- 
ence 1 are presented herein as a function of flow density (or mass ratio) 
for 12 wings with sweep angles from 0° to 52.5°, aspect ratios from 2.4 
7’b, taper ratios of 0.6 and 1.0, and center-of -gravity positions 
between 34 percent chord and 59 percent chord. The results of these 
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.palculations are examined to determine whether and under what conditions 
flutter at constant dynamic pressure is indicated. 

e 

In order to correlate data obtained under varying test conditions, 
the results of many recent flutter Investigations (refs. 6 to 10, for 
example) have been presented as the ratio of experimental flutter speed 
to that calculated by using two-dimensional Incompressible-flow aerody- 
namic parameters. As indicated in references 1 and 2, it has been pre- 
viously assumed that the resulting flutter-speed ratio was only slightly 


dependent on density. This supposition presupposes that density changes 
will affect experimental flutter speeds and those calculated from two- L 

dimensional incompressible-flow aerodynamics in approximately the same 4 

way. Reference 5 indicated that this supposition was not valid fqr low- 6 

aspect-ratio swept wings at transonic Mach numbers. In the present 4 

report the dependence of the flutter-speed ratio on flow density at sub- 
sonic and supersonic Mach numbers is examined by using the flutter cal- 
culation procedure of reference 1 for the 12 wings mentioned previously. 


This report also considers a slight modification to the method of 
reference 1 with respect to the location of the local aerodynamic center 
at subsonic speeds. Specifically, it is suggested that when calculated 
subsonic distributions of section lift-curve slope are used, the local 
aerodynamic centers be taken at the quarter chord of streamwise sections 
rather than at the quarter chord of sections normal to the elastic axis. 

The noncirculatory aerodynamic terms appearing in the flutter cal- 
culation procedures of references 1, 11, and 12 have been thought to have 
only a small effect on the calculated flutter speed. The importance of 
these terms is examined herein by comparing experimental flutter data and 
flutter calculations made by the complete method of reference 1 with 
flutter calculations for several wings made by omitting the noncircula- 
tory aerodynamic terms from the method of reference 1. 


SYMBOLS 


A aspect ratio of full wing including fuselage intercept 

a nondimensional distance from midchord to elastic axis 

measured perpendicular to elastic axis, positive rear- 
ward, fraction of semichord b 

ac nondimensional distance from leading edge to local aero- 

dynamic center (for steady flow) measured streamwise, 
fraction of streamwise chord 
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ac n 


ac n, l 
Aa c n 
b 
b r 

C 

C l 
C i 

6 oc, n 
C l 

L a,n, i 
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F 

G 

h 


^nr 


M 


M 


n 


Ma 
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nondimens Ional distance from midcbord to local aerodynamic 
center (for steady flow) measured perpendicular to elastic 
axis, positive rearward, fraction of semichord b 

values of ac n obtained from linearized aerodynamic theory 

change of ac n 

semi chord of wing measured perpendicular to elastic axis 

semichord of wing measured perpendicular to elastic axis at 
spanwise reference station tj = 0.75 

complex circulation function, F + iG 

local lift-curve slope for a streamwise section in steady 
flow 

local lift -curve slope for a section perpendicular to 
elastic axis in steady flow 

values of C ? obtained from linearized aerodynamic 

nx,n 

theory 

multiplier used to obtain variations in magnitude (but not 
in distribution) of ^ 

circulation function which modifies in-phase load components 

circulation function which introduces out-of-phase load 
components 

local vertical translational displacement of wing at elastic 
axis 

reduced frequency based on spanwise reference station 
(q = O.75) and on velocity component normal to elastic 
axis, b xU>/v n 

Mach number 

Mach number component perpendicular to leading edge 

oscillatory moment about elastic axis per unit length of 
wing, positive leading edge up 
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m r 
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Vr 
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K r 


a c/4 

A ea 
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n 

0 

P 

CJ 
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0 ) 
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mass of wing per unit length along elastic axis evaluated 
at spanwise reference station (q = O.75) 

oscillatory lift per unit length of wing along elastic 
axis, positive downward 

flutter speed, measured parallel to free stream (experi- 
mental values or values calculated by method of ref. l) 

calculated reference flutter speed obtained by using 
c l aj n = 2rt ' ac n = - an <i C = Fj + iGj 


component of free-stream velocity perpendicular to elastic 
axis 


nondimens ional distance from elastic axis to local center 
of gravity measured perpendicular to elastic axis, posi- 
tive rearward, fraction of semichord b 


wing section mass ratio evaluated at spanwise reference 


station (q = O.75) > 




m-v 


sweep angle of wing quarter chord, positive for sweepback 
sweep angle of wing elastic axis, positive for sweepback 

taper ratio of full wing including fuselage intercept 

nondimens ional coordinate (either spanwise or along elastic 
axis) measured from wing root, fraction of exposed panel 
span or fraction of panel length along elastic axis 


local torsional displacement of wing measured about elastic 
axis 


air density 

local bending slope of elastic axis 

local rate of change of twist angle with distance along 
the elastic axis 

circular frequency of vibration at flutter 

circular frequency of first uncoupled torsional vibration 
mode of wing measured about elastic axis 
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circular frequency of uncoupled bending vibration mode of 
wing (subscripts 1 and 2 denote first and second bending 
modes) 

circular frequency of vibration at flutter calculated by 

using = 2 jt, ac n = - and C = Fj + iGj 

cl iH ^ 


Subscripts: 

C circulation functions for compressible flow obtained in 

reference 1 from loading coefficients for two-dimensional 
wings oscillating in compressible flow 

I circulation functions for incompressible flow obtained in 

reference 11 for two-dimensional wings oscillating in 
incompressible flow 


PROCEDURE FOLLOWED IN THE F HITTER ANALYSIS 


Wing Designation 

The three-digit system used to identify the wings with taper ratio 
of 0.6 is the same as that used in reference 1. The first digit in this 
system is the aspect ratio of the full wing to the nearest integer. The 
second and third digits give the quarter-chord sweep angle to the nearest 
degree. For example, wing 445 has an aspect ratio of 4, a sweep angle 
of 45°, and a full-wing taper ratio of 0.6. Since some of the wings dis- 
cussed in this report have identical plan forms but different center-of- 
gravlty positions (see, for example, ref. 7), a single letter is appended 
to the plan-form designation to signify a forward or rearward shift in 
center of gravity. For example, wing 445 has a center of gravity at 
approximately 46 percent chord, whereas the center of gravity of wing 445F 
is at about 34 percent chord, and that of wing 445R is at about 58 percent 
chord. Wing 400 has a center of gravity at approximately 45 percent chord, 
but wing 400R has a center of gravity at about 59 percent chord. 

For the wings with taper ratio of 1.0, the same system is used, 
except that a fourth digit 1 is added to distinguish the taper ratio. 

For example, wing 4451 has a full-wing aspect ratio of 4, a sweep angle 
of 45°, and a taper ratio of 1.0. 

All of the wing plan forms treated in this report are shown in 
figure 1. 
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Determination of the Local Aerodynamic Center 
for Subsonic Speeds 

The flutter calculation procedure of reference 1 utilizes geometric, 
structural, and aerodynamic quantities associated with wing sections nor- 
mal to the elastic axis. In particular, for computational convenience, 
values of the section lift-curve slope C 7 and the local aerodynamic 

center ac, which are usually available for streamwise sections, are 

converted to C 7 and ac^ for sections normal to the elastic axis 

L a,n n 

by means of the geometrical relations 




L a,n cos 


(1) 


and 


ac n = 2ac \ cos A ea + 


tan A 


1 1 - A 


C A " A 1 + A 


(1 - 2 a) 


sin 


(1 + a) 


cos A„ + [tan A /. + i j- . ~ -- M sin A„„ 
ea \ c/4 A 1 + A/ ea 


I cos Ag a + a 


( 2 ) 


in the notation of the present report. 

For subsonic speeds below critical Mach number, it is reasonable to 
take the local aerodynamic centers at the local quarter-chord positions . 
However, the question immediately arises whether the aerodynamic center 
should be assumed to be at the quarter chord of a streamwise section 
|ac = or of a section normal to the elastic axis ^ac n = It may 

be noted that these two locations normally coincide for wings which are 
either unswept or untapered. The strip-theory flutter calculation method 
of reference 12 takes the aerodynamic center at the quarter chord of a 
section normal to the elastic axis. This location arose from the applica- 
tion of aerodynamic expressions for two-dimensional incompressible flow 
(from ref. 11) to strips normal to the elastic axis. In the present 
method, however, the local aerodynamic centers at subsonic speeds are 
taken at the quarter chord of streamwise sections for the following 
reasons: first, for subsonic speeds the spanwise distributions of C j 

are calculated by the lifting-line method of reference 13 as adapted for 
calculation in reference 14. The method of references 13 and 14 locates 
the lifting line and hence the aerodynamic center at the local streamwise 
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quarter chord ^ac = Aj. Second, the use of the streamwise quarter chord 

is consistent with the procedure employed for supersonic speeds; that is, 

C 7 and ac are determined for streamwise sections and then converted 

L <x 

t0 n and ac n "by means of equations (l) and ( 2 ). Thus the use of 

ac = i at subsonic speeds is consistent both with the method of calcu- 
lating at subsonic speeds and with the method for calculating 

and ac at supersonic speeds. 


It should be observed that for subsonic Mach numbers higher than 
critical it may become necessary to employ experimentally determined 
steady-flow aerodynamic parameters in order to account for the variation 
of aerodynamic center with Mach number. (See ref. 2.) 


For normally tapered (A < l) swept wings , the use of ac = — in 

in- 
equation ( 2 ) will yield ac n values of slightly greater magnitude (aero- 
dynamic center nearer the leading edge) than ac n = - A. The use of 

1 c - 

ac n = - — would thus result in slightly higher flutter speeds than the 
use of ac = A. As shown in reference 2, the use of incompressible-flow 

circulation functions (C = Fj + iGi) for nonzero subsonic Mach numbers 
results in flutter speeds slightly lower than those obtained by using 
C = ^(Fj + iGj), the form of circulation function appearing In the 


flutter calculation method of reference 1. For purposes of comparison 
some subsonic flutter calculations made by using ac n = -A and 

C = Fj + iGj are shown herein for the wings of figure 1 in addition to 

the calculations made by using ac = A 

4 

nitude of the factor Fq/Fj is indicated in figure 2. 


and C — jA^Fj + IGj^ • The mag— 


Systematic Variation of Aerodynamic Parameters and Density 

Variation of local aerodynamic center and local lift -curve slope . - 
As mentioned previously, the variations of n and ac n used herein 

were employed in two sets of subsonic and supersonic flutter calculations 
for three swept wings and four unswept wings (wings 445, 445F, 445R, 400, 
400R, 4001, and 7001) . In the first set of calculations 
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c l = c l 
L a,n L (x,n,l 


and 


ac n = ac n ,i + Aac n 


(3) 

(4) 


where /^ac n Is constant in any one flutter calculation but Is systemat- 
ically varied for each wing at each Mach number. This set of calculations 
thus indicates the effect on calculated flutter speed of arbitrary 
constant-percent -chord variations In local aerodynamic center. In the 
second set of calculations 


and 


= eC 


— CUT 

ct,n l 'tx,n,l 


(5) 


ac n = a c n ,l (6) 

where e is constant in any one flutter calculation but is systematically 
varied for each wing at each Mach number. This set of calculations thus 
Indicates the effect on calculated flutter speed of arbitrary variations 
in the magnitude (but not the distribution) of local lift-curve slope. 
Throughout these calculations flow density Is maintained constant for 
each wing. 


Variation of flow density .- As indicated previously, subsonic and 
supersonic flutter calculations have been made using various flow densi- 
ties for the 12 swept and unswept wings of figure 1. These calculations 
employed Cj ^ an( ^ ac n values obtained from linearized steady-flow 

aerodynamic theory. The density range covered for each wing generally 
encompasses the densities associated with all experimental flutter points 
given for the particular wing in references 6 to 10. 


Omission of noncirculatory aerodynamic terms . - In the flutter cal- 
culation procedure of reference 1 the section lift P and pitching 
moment Mq, are given by 


2 * * * • /,, . ' 

P = -npb h + v n 9 + v n d tan A ea - ba(0 + v n T tan 


- C l n n P v n bC Q 


Noncirculatory 

Circulatory 


( 7 ) 


and 


£r -P- fH 
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M a - -*pb^i + a 2 j(e + v n j tan Ag a ) 

+ jtpb v n ^h + v n o tan Ag a ^ + jtpb^a^h + v n o tan Agg J 
+ jtpb 2 v n 2 (e - abT tan Ag a j 


r Noncirculatory 


- 2jtpv n b 


I - ( a - “n)e 


Q 


Circulatory 


(8) 


where 


Q - h + v n 6 + v n a tan Ag a + b(-l 3 uJ 2 + a c n - aj(e + v n r tan A^) 

Since the noncirculatory terms of equations (7) and (8) essentially 
represent virtual-mass type effects, and since the method of reference 1 
is considered to be applicable only to cases involving low to moderate 
reduced frequencies, it has been thought that these noncirculatory terms 
have only a small effect on the calculated flutter speed. (See discussion 
of noncirculatory flow in refs. 1 , 15, and l6.) In order to investigate ‘ 
the importance of the noncirculatory terms, several subsonic and super- 
f° n * c ^ Utt f r calculatlons have been made for five swept wings (wings 445 
5 ^ 5 j ^52 ; and 430 ) by omitting the noncirculatory terms from the 

method of reference 1. These calculations are compared with calculations 
made by using the complete expressions for P and M a (eqs. (7) and ( 8 )) 
and with experimental flutter data. 

Vibration modes employed in the flutter calculations .- Since calcu- 
lated flutter characteristics appear not to be very sensitive to slight 
changes in mode shapes, and since the wings of this study are not highly 
tapered, flutter calculations for all of the wings presented herein were 
made by using the uncoupled mode shapes of a uniform cantilever beam. 

For all wings, the first torsion mode and the first and second bending 
modes were used. However, for wing 445 , subsonic and supersonic flutter 
calculations have also been made by using only the first torsion mode 
and the first bending mode. Comparison of the results of these two-mode 
calculations with the results of the three-mode calculations and with the 
experimental flutter data should give some indication of the importance 
of including higher modes in the analysis for both subsonic and super- 
sonic Mach numbers, at least for this particular wing. 
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RESULTS AND DISCUSSION 


Presentation of Results 


Subsonic and supersonic flutter calculations have been made by the 
method of reference 1 for 12 swept and unswept wings in order to evaluate 
the effects on calculated flutter characteristics of variations in aero- 
dynamic parameters and flow density. The plan forms of these 12 wings 
are shown in figure 1. Spanwise distributions of c i a ^ n = C l ajnj i an(i 

ac n = ac n z as obtained from linearized steady-flow aerodynamic theory 

for each wing are shown in figures 27 to 35 reference 1. These values 

0 f q 7 and ac 7 formed the aerodynamic basis for the present 

l a,n, l n U 

study. (See eqs. (3) to (6).) 

Structural data as well as experimental flutter data for all of these 
wings were obtained from references 1 and 6 to 10. The experimental 
flutter points shown herein were obtained at various values of density Pj 
whereas, for a particular wing, all of the points calculated in this 
investigation were obtained at a constant value of p (unless otherwise 
indicated) which represented approximately an average of the experimental 
densities. For each experimental point, however, the normalizing Vp 
was calculated by using the appropriate experimental density. The refer- 
ence flutter speeds used in references 6, 9> 10 l° r wings 245, 430, 

400 4001, and 7001 were calculated by employing only two degrees of 
freedom (first bending mode and first torsion mode). Since calculations 
employing three degrees of freedom (first and second bending modes and 
first torsion mode), as used herein, yield values of Vr which are 
slightly different from the two-degree-of-freedom values, the experimental 
y/Vr values for the wings listed above have been multiplied by the ratio 

Vr (for two degrees of freedom^ so that both calculated and experi- 
V R (for three degrees of freedom) 

mental flutter-speed, ratios as presented herein are normalized ty V R 
for three degrees of freedom. 

Variation of aerodynamic p arameters . - The effects on calculated 
flutter characteristics (flutter speed and frequency) of constant- 
percent-chord shifts in local aerodynamic center (eqs. (3) and (4)) are 
shown in figures 3 to 16 and are compared with experimental flutter data 
in figures 17 to 24. Figures 3 to lo also show the effects on calculated 
flutter characteristics of changes in the magnitude (but not in the dis- 
tribution) of local lift-curve slope (eqs. (5) and (6)). 

An empirical flutter-speed formula given in reference 17 indicates 
that, if the ratio of first bending -mode frequency to first torsion-mode 


-F” ctn -p- tr 1 
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frequency is small, then for a given wing at a given density the following 
proportionality applies 



(9) 


where C T is the wing lift-curve slope and - (f) „ is the chord ” 

wise distance between the section center of gravity and the section aero- 
dynamic center. If (J) - (^) aQ is evaluated at the reference station 

q = 0.75, then the proportionality relation ( 9 ) ma y l* e written in the 
notation of the present report as 

I a I i (10) 

Vr ^ + Xq, - ac n - 

From the relation (10) it appears that the principal effects of variations 
in lift -curve slope factor e might he accounted for by plotting je — 

instead of — . The relation (lO) also indicates that the principal 

Vr 

effects of constant -percent -chord variations in aerodynamic center Aac n 
might be accounted for by plotting ^a + % - ac n - Aac n )^_ 0 ^ ^ 
instead of In order to examine the possibilities of such correla- 

Vr 

tions the curves of figures 3 to 9 have teen replotted in figures 25 
to 31 as fe against e and - ac n - Aac n ^ =Q ^^ against 

Aac n . 


Figures 32 to 47 contain both experimental flutter data and subsonic 

flutter characteristics calculated by using ac = i- and C = ^(^I + d ®l) 

and by using ac n = - 1 and C = Fj + iGp. Figures 32, 3^, 36, 37, 38, 

40 42 44 45, and 46 also contain some flutter results obtained by 
omitting the noncirculatory aerodynamic terms from the flutter calcula- 
tion method of reference 1. 

Variation of flow density .- The variation of flutter-speed coeffi,- 
cient V/b^ with flow density p and mass ratio ic r is shown in 
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figures 48 to 58 . In these figures V/b^ (for constant Mach number) 
and p are plotted on logarithmic scales to facilitate determination of 
whether flutter at constant dynamic pressure is Indicated. The influence 
of density changes on flutter— speed ratio and flutter— frequency 

ratio o)/t% is illustrated in figures 59 to 80. From these calculations 
the amount of scatter to be expected in the experimental flutter data of 
references 6 to 10 based solely on density changes has been determined. 
The results are compared with the actual experimental scatter in fig- 
ures 8 l to 104. 

Finally, in figures 105 and 106 flutter characteristics calculated 
by using two degrees of freedom (first torsion mode and first bending 
mode) are compared with flutter characteristics calculated by using 
three degrees of freedom (first torsion mode and first and second bending 
mode) and with experimental flutter data for wing 445 . 

An Index to the above-mentioned figures is provided In table I. 


Influence on Calculated Flutter Characteristics 
of Variations in Aerodynamic Parameters 

Local aerodynamic center .- As expected, figures 3 to 16 illustrate 
the following qualitative results: The variations of calculated flutter 

speed with systematic constant -percent-chord variations in local aero- 
dynamic center indicate that for all wings and Mach numbers calculated 
the flutter speed increases as the aerodynamic center moves rearward. 

As Mach number increases through the transonic range, causing the aero- 
dynamic center to move rearward, the flutter speed becomes increasingly 
sensitive to changes in aerodynamic center. Further, comparison of 
results for wings 445, 445F, and 445R (figs. 3 to 5 ) and. for wings 400 
and 400R (figs. 6 and 7 ) indicate that the farther rearward the local 
center of gravity, and hence the greater the distance between center of 
gravity and aerodynamic center, the less the sensitivity to aerodynamic- 
center changes. This statement applies both to the systematic variations 
of aerodynamic center Aac n and to changes in ac n resulting from 
changes in Mach number . It may be observed from figures 6 and 9 that in 
comparison with results for the other wings, the calculated flutter speeds 
for wings 400 and 7001 are extremely sensitive to aerodynamic-center 
changes at supersonic Mach numbers. This great sensitivity arises because 
for these two wings at supersonic speeds the local aerodynamic centers 
are very close to the local centers of gravity. 

The calculated flutter frequencies shown in figures 10 to 16 show 
that with one exception the flutter frequencies increase as aerodynamic.^ 
center moves rearward and that as Mach number increases, the frequency 
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as well as the flutter speed becomes increasingly sensitive to aerodynamic 
center changes. The exception to this statement is wing 7001 at super- 
sonic speeds. Figure 16 shows that flutter frequencies for wing 7^01 
at supersonic speeds actually decrease slightly as aerodynamic center is 
moved rearward. 

The effects of systematic changes of aerodynamic center on the vari- 
ation of flutter speed and frequency with Mach number are shown in the 
crossplots in figures 17 to 24. Figures 17, 18, 21, and 22 show again 
that the rearward center of gravity of wing 400R results in variations 
(with ac n ) of flutter speed and frequency which are much smaller than 
those associated with wing 400. Inasmuch as linearized aerodynamic 
theory characteristically predicts aerodynamic centers which are farther 
rearward than those obtained experimentally, it would appear that flutter 
calculations using linearized-theory aerodynamics should be more reliable 
for wing 400R than for wing 400 especially at supersonic speeds. The 
solid curves of figures 17, 18, 21, and 22 show this to be true. Fig- 
ure 20 shows that for wing 7001 good agreement between calculated and 
experimental supersonic flutter speeds Is obtained with local aerodynamic 
center only one percent chord forward of the linearized-theory position 
(Aac n = -0.02) . 

Quantitatively, for wings 400 and 400R results obtained by using 
linearized-theory aerodynamic parameters may be compared with results 
obtained by using Cj and ac n from flight tests of the X-1E air- 
plane since the wing geometry of the X-1E was closely similar to wings 400 
and 400R. Figure 5 of reference 2 shows that over the middle and outboard 
portions of the wing panel C ia n at M = 1.41 for the X-1E airplane is 

very close to Cq at M = ^2* obtained from linear theory. However, 

the local aerodynamic centers for the X-1E are about 6.3 percent chord 

forward of those given by linear theory. Use of C7 and ac n for 

L ct,n 11 

the X-1E in flutter calculations for wings 400 and 400R at M = (2 
results in a 51-percent reduction of flutter speed for wing 400 but only 
a 9-P erc e n t reduction for wing 400R compared with linear-theory results. 
(See figs. 6, 7; 17 > and- 18.) Similarly use of n and ac n for 

the X-1E results in a 54-percent reduction of calculated flutter frequency 
for wing 400 but only a 2-percent reduction for wing 400R. (See figs. 13, 
14, 21, and 22.) These comparisons give some indication of the magnitude 
of error which can occur when linear-theory aerodynamic parameters are 
used in flutter calculations for finite-thickness wings in a real fluid, 
particularly when the local aerodynamic centers are close to the local 
centers of gravity. 


In figures 32 to 47 subsonic flutter characteristics calculated by 
using ac n = - i and C = Fj + iGj are compared with results obtained 



1 6 


and C = — (Fj + iGj). As expected, the flutter speeds 


ty using ac = - 

and frequencies calculated by these two procedures generally differ by 
only a very small amount. For the swept, tapered wings the maximum 
difference between flutter speeds calculated by these two procedures 


occurs at M = 0 


where — = 1.0, and 
F I 


ac n = - i yields higher flutter 


speeds than ac = i. The largest difference (about 4 percent) occurs for 

wing 445F because the forward local centers of gravity of that wing are 
comparatively close to the local aerodynamic centers, and that condition 
yields increased sensitivity to ac n changes. (Compare figs. 3, 4, 
and 5.) As Mach number increases from 0, Fq/Fj decreases from 1.0, 

and subsonic flutter speeds calculated by the two above-mentioned proce- 
dures come closer together. For all of the swept, tapered wings except 
wings 245 and 445F the difference is insignificant at M = 0.75* For 
the three untapered wings (wings 4451, 4001, and 7001) substitution of 
ac = i into equation (2) yields ac n = -i so that the two above- 

4 ^ , 

mentioned procedures give identical results at M = 0. (See figs. 19, 

20, 23, 24, 39, and 47.) At higher subsonic Mach numbers the two proce- 
dures reflect only the difference between using C = + iGj) and 


C = Fj + iGj. As indicated in references 1 and 2 this difference should 
be small as long as the Mach number component normal to the leading edge 
is not large. This statement is supported here by figures 19, 20, 23, 

24, 39, and 47. For swept wing 4451 at M = 0.75 (fig- 39) the two 
procedures yield flutter speeds which differ by only 2 percent. For 
unswept wing 7001 (fig. 20), however, the component of Mach number normal 
to the leading edge is the same as stream Mach number, and at M = 0.75 
the two procedures yield flutter speeds which differ by about 8 percent. 
The corresponding differences in flutter frequency, however, are insig- 
nificant for all of the wings. 


Figures 25 to 31 show attempted correlations of the effects of ac n 

changes based on the relation given by equation (10) . Most of the curves 
of figures 25 to 31 indicate little variation of the quantity 

Qa + Xa - ac n - Aac n )^_ Q ^ ^ with changes in ac n . The exceptions 

to this statement are the curves for the highest Mach numbers and for the 
wings with the more forward local centers of gravity. Thus it appears 
that equations (9) and (10) should predict with reasonable accuracy the 
effects of ac n changes as long as the changes Aac n are not too large 
and as long as the local aerodynamic centers are not close to the local,, 
centers of gravity. 
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. ,, L ^ Ca , 1 V ' 1 j~ ft ~^ UrVe slo P e -- As expected (see eq. ( 10) ) , figures 3 to 9 
indicate that calculated flutter speed generally decreases as local lift- 
curve slope increases. The one exception to this statement is the curve 
shown m figure 9 for wing 7001 at M = 1.30. This curve shows an unex- 
pected increase of flutter speed as lift-curve slope increases. An exam- 
ination of the structural and aerodynamic quantities for wing 7001 
revealed that the increase of flutter speed with increasing lift-curve 
slope is associated with the location of local aerodynamic centers rear- 
ward of the elastic axis. For wing 7001 at M = I.30, this condition 
existed over the major portion of the wing and resulted in a change of 
sign for the circulatory component of the section pitching moment M^ 

(Se e eq (8).) For wing 7001 at M = I.I547O and for wing 400 at 
V"*!* 1 , ° and M = \[2, the condition existed over a smaller portion 
o the wing and resulted in comparatively small decreases of flutter 
speed with increasing lift-curve slope. (See figs. 6 and 9.) 


Figures 25 to 31 show attempted correlations of the effects of 


changes phased on the relation given by equation (lO). Most of the curves 
° f Vr show a tendency to rise with increasing e, indicating that 

calculated flutter speed generally does not decrease quite as rapidly as 
the inverse square root of lift-curve slope. The deviations previously 
described for wings 400 and 7 001 at supersonic speeds are clearly brought 
° ut in ^Sures 28 and 31. Figures 25 to 31 seem to indicate that equa- 
10ns (9) and (10) should predict with reasonable accuracy the effects 
C £a,n c h an g e s as long as the local aerodynamic centers are not close 

to the local centers of gravity and as long as the local aerodynamic 
centers are forward of the elastic axis. 


'a,n 


Influence on Calculated Flutter Characteristics 
of Variations in Flow Density 

In figures 48 to 58 the calculated flutter-speed coefficient — 
for 11 wings has been plotted logarithmically as a function of flow 


density and as a function of the mass ratio k 


itph 2 


rru 


On these plots 


flutter at constant dynamic pressure would be indicated by a straight 
line of slope - These figures show that for most of the wings calcu- 
lated the variation of flutter-speed coefficient — ¥ — with flow density 

br<%, J 

for nonzero Mach numbers has the same general character as the variation 

vp w 

° f bpi\x' Also > for a11 °f the swept wings the curves of figures 48 to 58 
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have slopes of magnitude significantly less than I— — J. However, for 

V R V 

unswept wings 400 and 7^04 the curves of and of — for 

M = 0.75 closely follow a slope of - indicating flutter at approxi- 
mately constant dynamic pressure. For unswept wings 400R and 4001 all 
of the curves shown indicate flutter at approximately constant dynamic 

pressure. The approximate slope of - i for f° r a -^ unswe P b 

wings supports the statement in reference 5 that for the V R calcula- 
tion, flutter at constant dynamic pressure is indicated only for unswept 
wings or for swept wings of high aspect ratio. Figures 54 and 58 show, 

however, that even though =r~- values for the unswept wings indicate 

b r°tt 

flutter at constant dynamic pressure, constant-dynamic-pressure flutter 
may not occur for supersonic Mach numbers, particularly when local aero- 
dynamic centers lie close to local centers of gravity as occurs for 
wings 400 and J001 . 

The influence of density changes on calculated flutter-speed ratio 
v/Vr and flutter-frequency ratio c 0 / 1 % is illustrated in figures 59 
to 80. In figures 59 to 69 the small variations of calculated subsonic 
Y^Yr with density changes seem to indicate that the presentation of sub- 
sonic flutter speeds in the form V/V R is probably acceptable if the 
density range covered is no greater than those shown. Use of the V^V R 

form for supersonic flutter speeds, however, may result in unacceptable 
scatter if the density range covered is large. 

Figures 70 to 80 generally show appreciable increases of flutter 
frequency with increasing density and thus indicate significant increases 
in the aerodynamic contribution to the apparent wing stiffness. The 
exceptions to this statement are wing 400 at M = {2 (fig. 76 ) and 
wing 7001 at M = 1.15470 and M = 1.30 (fig. 80) . These are again 
the conditions for which the local aerodynamic centers lie rearward of 
the elastic axis over a large portion of the wing. 

From figures 59 to 80 the amount of scatter to be expected in the 
experimental flutter data of references 6 to 10 based solely on density 
changes has been determined. Figures 8 l to 104 show comparisons of these 
estimates and the actual scatter of the data. These figures show that 
for some of the wings density changes may be responsible for an appreci- 
able portion of the scatter. It should also be noted that the experi- 
mental scatter shown in these figures tends to be systematic, and the 
experimental flutter-speed ratios tend to vary with density in the same 
manner as the calculated values, particularly at supersonic Mach numbers. 


CN 
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Omission of Noncirculatory Aerodynamic Terms 

As indicated earlier, some flutter calculations have been made with 
the noncirculatory components of section lift P (eq. (7)) and pitching 
moment (eq. (8)) omitted. The results of these calculations 

(figs. 32, 34, 36, 37, 38, 40, 42, 44, 45, and 46) indicate that the 
effects of the noncirculatory terms on subsonic and supersonic flutter 
speeds are small (of the order of 15 percent or less) but not negligible. 
The effects on flutter frequencies of omitting the noncirculatory terms 
are appreciably greater than the effects on the flutter speeds. Thus the 
importance of the noncirculatory terms is not insignificant, and they 
should not be deleted from subsequent flutter calculations. 


Effect of the Higher Modes on Calculated 
Flutter Characteristics 

As mentioned earlier, all of the flutter calculations previously 
presented were made by using the first torsion mode and the first and 
second bending modes of a uniform cantilever beam. In order to obtain 
a preliminary estimate of the importance in the flutter analysis of the 
second and higher modes, subsonic and supersonic flutter calculations 
have been made for wing 445 hy using only the first torsion mode and 
the first bending mode. These calculations are compared with the three- 
mode results in figures 105 and 106. It is evident that the second 
bending mode contributes significantly to the flutter-speed ratio 
only at the higher supersonic Mach numbers shown. The effect of the 
second bending mode on the actual speed, however, is indicated by the 
fact that 


V R ( two modes) 
Vr (three modes) 


1.038 


so that at M = O.75, 

V (two modes) 
V (three modes) 


1.051 


but at M = 1.75, 

V (two modes) _ ^ ^g, 

V (three modes) 

Flutter frequencies calculated with two modes and with three modes differ 
by only a small amount through the Mach number range calculated. Thus, 


S 
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as might be expected, when flutter frequencies increase (as Mach number 
increases, for example) the higher modes become more important in the 
flutter calculation. 


CONCLUSIONS 


Subsonic and supersonic flutter calculations have been made by the 
method of NACA Research Memorandum L57L10 for 12 swept and unswept wings 
in order to evaluate the effects on calculated flutter characteristics 
of variations in aerodynamic parameters and flow density. For all of 
these wings the ratio of first bending-mode frequency to first torsion- 
mode frequency was of the order of 0.3 or less. With section lift-curve 
slopes and local aerodynamic centers obtained from linearized aerodynamic 
theory used as the basis, systematic constant-percent-chord shifts of 
local aerodynamic center and systematic changes in the magnitude (but 
not in the distribution) of section lift-curve slope indicate the 
following: 

1. If the local aerodynamic centers are sufficiently forward of the 
local centers of gravity, the calculated flutter speed appears to vary 
approximately as the inverse square root of the distance between local 
aerodynamic center and local center of gravity . 

2. When local aerodynamic centers are close to the local centers of 
gravity, calculated flutter speed is extremely sensitive to small changes 
in aerodynamic-center position, and the above empirical rule is no longer 
valid. When this latter condition exists, flutter calculations based on 
linear -theory aerodynamics are likely to be considerably in error because 
linearized theory characteristically predicts aerodynamic centers that 
are farther rearward than those observed on finite-thickness wings in a 
real fluid. 

3. If the local aerodynamic centers are forward of the elastic axis 
and are not close to the local centers of gravity, the calculated flutter 
speed appears to vary approximately as the inverse square root of lift- 
curve slope . 

4. If local aerodynamic center moves rearward of the elastic axis, 
a change of sign occurs in the circulatory component of section pitching 
moment, and calculated flutter speeds may actually increase with increasing 
lift-curve slope. 

Some variations of flow density employed in flutter calculations 
based on steady-flow linearized-theory aerodynamics indicate the following: 


4r- Q> 
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5 . Flutter at constant dynamic pressure (flutter speed varying as 
the - i power of density) was indicated only for unswept wings for which 

the local aerodynamic centers were sufficiently forward of the local 
centers of gravity. ... 

6 . For unswept wings with aerodynamic centers close to the centers 
of gravity and for all of the swept wings, calculated flutter speed ^ 
varied approximately as a power of density of magnitude smaller than - — 

7 . Presentation of flutter speeds in the form of a flutter-speed 
ratio, 'flutter speed divided by flutter speed calculated from two- 
dimensional incompressible -flow aerodynamics, is probably acceptable 
for most wings at subsonic Mach numbers. For supersonic Mach numbers, 
however the increased variation of this ratio with density changes may 
introduce unacceptable scatter if the density range covered is large. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , August 25, 1959 - 
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Figure 2.- Variation with reduced frequency of the factor Fq| Fj 

which accounts for the effect of compressibility in the circulation 

function C = — (Ft + iG-r'l . 
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Calculation using Ci a n a 
ac n obtained from flight 
test of the X-1E airplane 
atM = 1.1*1 (See ref. 2) 
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Figure 6.~ Variation of flutter speed with lift factor e and with aerodynamic-center shift 
A&c n for wing 400, V R is calculated for e = 1 and Aac n =0. p = 0.002378 slug/cu ft 
V R = 976.5 ft/sec. 












for wing 445- P = O.OO38OO slug/cu ft; a) a = 2,192 radians/; 



Figure 11.- Variation of flutter frequency with lift factor e and with aerodynamic-center shift 
Aac n for wing 445F. p = 0.003000 slug/cu ft; = 1,144 radians/sec. 






Figure 12.- Variation of flutter frequency with lift factor e and with aerodynamic -center shift 
Aac n for wing 445R. P = 0.002378 slug/cu ft; o> a = 2,306 radians/sec. 
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.gure 14.- Variation of flutter frequency with lift factor e and with aerodynamic -center shift 
£ac n for wing 400R. p = 0.003100 slug/cu ft; o> a = 1,982 radians/sec. 




Figure 15.- Variation of flutter frequency with lift factor e and with aerodynamic -center shift 
Aac n for wing 4001. p = 0.002378 slug/cu ft; c% = 2,048 radians/sec. 
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Figure 16.- Variation of flutter frequency with lift factor e and with aerodynamic -center shift 
Aac n for wing 7001. p = O.OO 55 OO slug/cu ft; = 2,271 radians/sec. 









Flutter calculation 

Aac 



Figure 18.- Effect of ac n changes on the variation of flutter speed with Mach number for 
wing hOOR. For calculated points p = 0.003100 slug/cu ft and V R = 852.5 ft/sec. 





Figure 20.- Effect of ac n changes on the variation of flutter 
speed with Mach number for wing 7001. For calculated points 
p = 0.005500 slug/cu ft and Vp = 844.8 ft/sec. 






Flutter experiment 
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Figure 21.- Effect of ac n changes on the variation of flutter frequency with Mach number for 
wing 400. For calculated points p = 0.002378 slug/cu ft and cu a = 2,463 radians/sec. 
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Figure 23.- Effect of ac n changes on the variation of flutter frequency with Mach number for 
wing 4001. For calculated points p = 0.002378 slug/cu ft and = 2,048 radians/sec. 
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igure 24.- Effect of ac n changes on the variation of flutter frequency with Mach number for 
wing 7001. For calculated points p = 0.005500 slug/cu ft and cjo a = 2,271 radians/sec-. 















with aerodynamic-center shift Aac n for wing 4-00. Vr is calculated for e = 1 and 
Aac n =0. p = 0.002378 slug/cu ft; Vp = 976.5 ft/sec. 


















subsonic Mach numbers 



Figure 32.- Variation of flutter speed with Mach number for wing h-h-5 . For calculated points 

p = O.OO 38 OO slug/cu ft and V R = 735.0 ft/sec. 







Flutter experiment 



0.002378 slug/cu ft and V R = 928.7 ft/; 
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Figure 35 .- Variation of flutter Bpeed with Mach number for 
wing 2^5. For calculated points p = 0.003900 slug/cu ft 
and Vp = 65 O .4 ft/sec. 



Flutter experiment 



0.005500 slug/ cu ft and V R = 901.1 ft/; 




Calculation usiryz; ac 



Figure 37 •- Variation of flutter speed with Mach number for wing 1+52. For calculated points 

p = 0.002700 slug/cu ft and Vp = 988.7 ft/sec. 
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Figure 38 .- Variation of flutter speed with Mach number for wing 430. For calculated points 

P = 0.003700 slug/cu ft and V R = 694.7 ft/sec. 




Calculation using 



Figure 59 *“ Variation of* flutter speed with Mach number for wing 4451 * For calculated points 

p = 0.003200 slug/ cu ft and Vr = 944.7 ft/sec. 





Calculation using ac * l/k 



0.003800 slug/cu ft and = 2,192 radians/; 
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Figure 41.- Variation of flutter frequency with Mach number for wing 445F. For calculated points 

p = 0.003000 slug/cu ft and u) a = 1,144 radians/sec. 














Platter experiment 



Figure 44.- Variation of flutter frequency with Mach number for wing 645. For calculated points 

p = 0.003500 slug/cu ft and oo a = 3 >173 radians/sec. 
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Figure 45 .- Variation of flutter frequency with Mach number for wing 452. For calculated points 

p = 0.002700 slug/cu ft and = 2,500 radians/sec. 




Figure 46.- Variation of flutter frequency with Mach number for wing 450. For calculated points 

p - 0,005700 slug/cu ft and = 2,158 radians/sec. 






Flutter experiment 



Figure 47-- Variation of flutter frequency with Mach number for wing 4451- For calculated points 

p = 0.005200 slug/cu ft and = 2,552 radians /sec. 
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Figure 48.- Variation of flutter-speed coefficient with density for 

wing 445. 
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Figure 49- - Variation of flutter- speed coefficient with density for 

wing 445F. 
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Figure 52.- Variation of flutter- speed coefficient with density for 

wing 1 + 52 . 














Figure 57.- Variation of flutter- speed coefficient with density for 

wing 1+001. 
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figure 59-- Variation of flutter-speed ratio with flow density for 

wing 445. 
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Figure 6l.- Variation of flutter-speed ratio with flow density for 

wing 44 5R. 
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Figure 62.- Variation of flutter-speed ratio with flow density for 

wing 645* 
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Figure 65 .- Variation of flutter-speed ratio with flow density for 

wing 400. 
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Figure 66.- Variation of flutter-speed ratio with flow density for 

wing 400R. 
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Figure 67 .- Variation of flutter-speed ratio with flow density for 

wing 4451. 
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Figure 69 .- Variation of flutter-speed ratio with flow density for 

wing 7001* 
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Figure 70 .- Variation of flutter frequency with flow density for 
wing 445- c% = 2,192 radians/sec. 
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Figure 72 .- Variation of flutter frequency with flow density for 
wing = 2,306 radians /sec. 
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Figura 74 • — Variation of 1 fluttar fracjuancy with fiow dansity for 
wing 452. ca a = 2,300 radian8/sec. 
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Figure 77.- Variation of flutter frequency with flow density for 
wing 400R. = 1,982 radians/sec. 
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Figure 78*- Variation of flutter frequency with flow density for 
wing 1+451. <% = 2,552 radians/sec. 
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Figure 79-- Variation of flutter frequency with flow density for 
wing 4001. = 2,048 radians/sec. 
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Figure 80.- Variation of flutter frequency with flow density for 
wing 7001. cDq, = 2,271 radians/sec. 
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Flutter experiment (p varies) 
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Figure 8l.- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing V and V-^ are both calculated for the density indicated. 



Flatter experiment (p varies) 



Figure . — Effect, of density changes on the variation of flutter— speed ratio with Mach number 
for wing 445F . V and V R are both calculated for the density indicated. 




Flutter calculation 



Figure 83 .- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing 44-5R. V and Vjj are both calculated for the density indicated. 








Figure 84.- Effect of density changes on the variation of flutter-speed 
ratio with Mach number for wing 245- V and are both calculated 

for the density indicated. 



Figure 85 .- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing 645- V and Vr are both calculated for the density indicated. 






Flutter experiment ,(p varies) 
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Flutter experiment (p varies) 
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Figure 87 .- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing k30. V and Vjj are both calculated for the density indicated. 
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Figure 88.- Effect of density changes on the variation of flutter-speed 
ratio with Mach number for wing 400. V and Vp are both calculated 
for the density indicated. 


Flutter experiment (p varies) 
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Figure 89 .- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing 100R. V and Vp are both calculated for the density indicated. 



Flutter experiment (p varies) 
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Figure 90.- Effect of density changes on the variation of flutter-speed ratio with Mach number 
for wing 4451. V and V R are both calculated for the density indicated. 
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•Figure 91 •- Effect of density changes on the variation of flutter-speed 
ratio with Mach number for wing 4001. V and Vp are both calcu- 
lated for the density indicated. 
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Figure 92.- Effect of density changes on the variation of flutter-speed 
ratio with Mach number for wing J001. V and V R are both calcu- 
lated for the density indicated. 










Flatter calculation 



Figure 9^-“ Effect of density changes on the variation of flutter frequency with Mach number for 

wing 1^5F. = 1,144 radians/sec. 
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Figure 96 .- Effect of density changes on the variation of flutter fre- 
quency with Mach number for wing 245- = 1*665 radians/sec. 
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Flutter calculation 



Figure 98*** Effect of density changes on the variation of flutter frequency with Mach number for 

wing 452. <% a 2,500 radians/sec. 



Flutter experiment (p varies) 



Figure 99* - Effect of density changes on the variation of flutter frequency with Mach number 

wing 430. = 2,158 radians/sec. 



Flutter experiment (p varies) 
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Figure 100.- Effect of density changes on the variation of flutter frequency with Mach number 

for wing 400. cn a = 2,463 radians/sec. 





Flutter experiment (p varies) 



Figure 101.- Effect of density changes on the variation of flutter frequency with Mach number 

for wing 400R. = 1,982 radians/sec. 










Flutter experiment (p varies) 



Figure 103.- Effect of density changes on the variation of flutter frequency with Mach number 

for wing 4001. o) a = 2,048 radians/sec. 




Flutter experiment (p varies) 



Figure 104*- Effect of density changes on the variation of flutter frequency with Mach number 

for wing 7001. c% = 2,271 radians/sec. 









